Absence of magnetic order for the spin-half Heisenberg antiferromagnet on the star 

lattice 
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We study the ground-state properties of the spin-half Heisenberg antiferromagnet on the two- 
dimensional star lattice by spin-wave theory, exact diagonalization and a variational mean-field 
approach. We find evidence that the star lattice is (besides the kagome lattice) a second candi- 
date among the 11 uniform Archimedean lattices where quantum fluctuations in combination with 
frustration lead to a quantum paramagnetic ground state. Although the classical ground state of 
the Heisenberg antiferromagnet on the star lattice exhibits a huge non-trivial degeneracy like on 
the kagome lattice, its quantum ground state is most likely dimerized with a gap to all excitations. 
Finally, we find several candidates for plateaux in the magnetization curve as well as a macroscopic 
magnetization jump to saturation due to independent localized magnon states. 

PACS numbers: 75.10.Jm; 75.45.-l-j; 75.60.Ej; 75.50.Eg 



I. INTRODUCTION 

The spin-half two-dimensional (2D) Heisenberg anti- 
ferromagnet (HAFM) has attracted much attention in 
recent times. In particular, the recent progress in syn- 
thesizing novel quasi-2D magnetic materials exhibiting 
exciting quantum effects has stimulated much theoreti- 
cal work. We mention for example the spin-gap behav- 
ior in CaV405^ and in SrCu2 (603)2,'^ the spin fraction- 
alization in Cs2CuCl4'^ or the plateau structure in the 
magnetization process of frustrated quasi-2D magnetic 
materials like SrCu2(B03)2^^ or Cs2CuBr4.'^ 

While the ground state (GS) of the one-dimensional 
quantum HAFM does not possess Neel long-range or- 
der, for the spin- half HAFM on 2D lattices the com- 
petition between quantum fluctuations and interactions 
seems to be well balanced and magnetically ordered and 
disordered GS phases may appear. A fine tuning of this 
competition may lead to zero temperature transitions be- 
tween semi-classical and quantum phases. The proto- 
types of 2D arrangements of spins are the 11 uniform 
Archimedean lattices (tilings), see, e.g. Refs. 5 6 7. They 
present an ideal possibility for a systematic study of the 
interplay of lattice geometry and magnetic interactions 
in 2D quantum spin systems. 

The HAFM on the widely known square, honeycomb, 
triangular and kagome lattices has been studied in nu- 
merous papers over the last decade. While for the square, 
honeycomb and triangular lattices the existence of semi- 
classical r nagn etic order seems to be well-established (see 
e.g. Refs. '7?'?) the spin-half HAFM on the kagome lat- 
tice is a candidate for a magnetic system with a quantum 
paramagnetic GS (see the reviews mlalglllOl and references 
therein) . Other less known Archimedean lattices like the 
maple-leaf latticcfii the square-hexagonal-dodecagonal 
latticeiSiH and the trellis latticei^ji^ exhibit most likely 
semi-classical magnetic GS order. 




FIG. 1: The star lattice with iV = 42 sites. 



In this paper we present another candidate for a quan- 
tum paramagnetic GS among the Archimedean lattices, 
namely the so-called star lattice, featured by low coor- 
dination number z = ?> and strong frustration due to a 
triangular arrangement of bonds (see Fig. ^ . 



II. MODEL 

The geometric unit cell of the star lattice contains six 
sites and the underlying Bravais lattice is a triangular 
one (see Figs. ^ and |2l. For this lattice we consider the 
spin-half HAFM in a magnetic field h 
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where the sum runs over pairs of neighboring sites (ij) 
and S^ — J2i^i- The star lattice carries topologically 
inequivalent nearest-neighbor (NN) bonds Jd (dimer 
bonds, solid lines in Fig. QJ and Jt (triangular bonds. 
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FIG. 2: (Color online) Two variants of the GS of the classical 
HAFM on the star lattice: the \/3xi/3 state (left) and the 
q=0 state (right). The dotted ellipses show further degrees 
of freedom of the highly degenerate classical GS. Different 
shades of the triangles symbolize different chiralities. 



dashed lines in Fig. ^ see also Fig. ^ . For the uniform 
lattice these bonds are of equal strength Jd ^ Jt = J ■ 



III. SEMI-CLASSICAL GROUND STATE 

In the classical GS for ft, = the two non-equivalent 
NN bonds of the star lattice carry different NN spin cor- 
relations: The angle between neighboring spins on dimer 
bonds Jd is tt, whereas the angle on triangular bonds Jt 
is 27r/3. Its energy per bond is e'^^^^ = —1/6. The classi- 
cal GS for the star lattice has a great similarity to that of 
the kagome lattice. It also exhibits a non-trivial infinite 
degeneracy. Moreover, there are also the two variants of 
the classical GS, namely the so-called \/3x \/3 and q — 
states (see Fig.|2l, often used for discussing possible or- 
der in the kagome lattice. Hence these two particular 
planar states can also be considered as variants of possi- 
ble GS ordering for the star lattice. In the following we 
discuss the influence of quantum fluctuations on the GS 
properties. 

First, we perform a linear spin- wave theory (LSWT) 
starting from the planar classical GSs. The LSWT for the 
star lattice is more ambitious than for the kagome lattice, 
since we have to consider six types of magnons. As in 
the kagome caseiStiLi^, the spin- wave spectra are equiva- 
lent for all coplanar configurations satisfying the classical 
GS constraint. We obtain six spin-wave branches. Two 
dispersionless modes are found, namely wiq — 0, Ci'2q = 
Jsv3. Thus also a fiat zero-mode appears as it is ob- 
served for the kagome case. In addition there are two 
acoustical and two optical branches. There is no order- 
by-disorder selection among the coplanar classical GSs 
due to the equivalence of the spin- wave branches obtained 
from LSWT, exactly like for the kagome latticeii2iii 
The GS energy per bond for s — 1/2 in LSWT is 
eo/J = —0.296759. Due to the fiat zero mode the in- 



tegral for the sublattice magnetization divergesiS. which 
might be understood as another hint for the absence of 
the classical order. Although on the semi-classical LSWT 
level both the kagome and the star lattice exhibit almost 
identical properties, the situation might be changed tak- 
ing into account the quantum fluctuations more properly. 



IV. LANCZOS EXACT DIAGONALIZATION 

We consider now the extreme quantum limit s = 1/2 
by direct numerical calculation of the GS and the low- 
lying excitations aX h = for finite lattices of A^ = 
18,24,30,36,42 sites. For each size N we have cho- 
sen only lattices having good geo metric properties us- 
ing the criteria given in Ref. ^^ The largest lattice 
{N — 42) is shown in Fig. ^ and required a Lanczos 
diagonalization in dimension 801 258 898 for the compu- 
tation of GS properties. The GS for all these lattices is 
a singlet and has an energy per bond eg = —0.312479 
{N = 18); -0.311342 (24); -0.310808 (30); -0.310348 
(36A); -0.310657 (36B); -0.309918 (42). The first exci- 
tation is a triplet and has a gap A = 0.578710 {N = 18); 
0.531822 (24); 0.498564 (30); 0.480343 (36A); 0.483112 
(36B) (no result available for TV = 42). 

We present in Tabled the spin-spin correlation for the 
largest finite lattice considered. Note that the two non- 
equivalent NN correlation functions differ drastically, we 
have (5'oS'lg) ~ 3.5(S'oS'f) indicating a tendency to form 
local singlets on the dimer bonds. 

Let us compare the spin correlations with those for the 
HAFM on triangular and kagome lattices. We consider 
the strongest correlations as a measure for magnetic order 
and present in Fig.|3|the maximal absolute correlations 
I {SqS^) |,„ax for a certain separation R = |Ro — Ri| versus 
R. As expected we have very rapidly decaying correla- 
tions for the disordered kagome case, whereas the corre- 
lations for the Neel ordered triangular lattice are much 
stronger for larger distances and show a kind of satura- 
tion for larger R. Although the correlations for the star 
lattice are slightly larger than those of the kagome lattice 
they are significantly weaker for separations i? > 3 than 



TABLE I: All non-equivalent GS spin-spin correlations 
{SqS!) = |{SoS») for the HAFM on the star lattice with 
A'' = 42 sites. In addition to the site index i we give the sep- 
aration R — |Ro — Ri| between sites and i in units of NN 
separation. 



i (R) 


1 (1) 

-0.05643 


3 (1.932) 
0.03559 


4 (2.909) 
-0.01058 


i (R) 


6 (3.732) 
-0.00451 


7 (3.346) 
0.01066 


9 (5.278) 
0.00439 


i (R) 


10 (4.625) 
-0.01173 


26 (2.732) 
-0.03875 


28 (1) 
-0.19707 



0.15 



a 
A_ 

N ' 

C/3 



C/3 

V 



0.05 



star N= 
kagome N= 
triangular H-- 



42 

:36 
:36 




5.5 




FIG. 3: Maximal spin-spin correlation KSoSf )|max versus 
separation R — |Ro — Ri| for the star lattice (A'^ = 42), the 
kagome {N = 36) and the triangular lattice (A^ = 36) (the 
lines are guides for the eyes) . The data for the kagome lattice 
coincide with those from Ref. 12(1 



those of the triangular lattice. The large NN correlation 
for the star lattice corresponds to a NN dimer bond. 

Next we consider the low-lying spectrum of the star 
lattice (see Fig. 0J, following the lines of the discus- 
sion of the spectrum for the triangular^-'^ and the kagome 
lattice. ^^'^^ It is obvious that the lowest states E^ainiS) 
are not well described by the effective low-energy Hamil- 
tonian Hcs ^ Eq + S'^/2Nxo of a semi-classically or- 
dered system: (i) The £^min versus S{S + 1) curve devi- 
ates significantly from a straight line (cf. the dashed line 
in Fig.0}. (ii) We do not see well separated lowest states 
in the different sectors of S (so-called quasi degenerate 
joint stated) which could collapse onto a Neel-like state 
in the thermodynamic limit, (iii) The symmetries of the 
lowest states in each sector of S cannot be attributed 
to the classical -s/SxVS or q=0 GSs in general. These 
features are similar to the kagome lattice. ^^-^"^ However, 
there is one striking difference. In contrast to the kagome 
lattice we do not have non-magnetic singlets filling the 
singlet-triplet gap (spin gap) commonly interpreted as a 
remnant of the huge classical GS degeneracy. Rather we 
have a particularly large singlet-singlet gap. This basic 
difference to the kagome lattice can be attributed to the 
special property of the quantum GS of the star lattice 
to have strongly enhanced antiferromagnetic correlations 
on the Jd bonds. As a consequence the quantum GS of 
the star lattice has an exceptionally low GS energy eg 
(see Table 111)1 and is well separated from all excitations. 

For finite systems the order parameter is based on the 
spin-spin correlation functions. For systems with well- 
defined semi-classical long-range order usually the square 
of the staggered magnetization is used. However, this 
definition of the order parameter is problematic in the 
present situation: due to the huge non-trivial degeneracy 
of the classical GS it remains unclear which type of order- 
ing might be favored in the quantum system. Therefore 



FIG. 4: Low-energy spectrum for the HAFM on the star lat- 
tice {N — 36B) (the inset shows the k points in the Brillouin 
zone). 



we use a definition of an order parameter 

1/2 



m 






•piSj) 



(2) 



which is independent of any assumption on classical 
order,- For bipartite systems like the square lattice this 
definition is identical to the staggered magnetization fh 
and for the HAFM on the triangular lattice (?tt,+)^ is by 
1/3 larger than the usual definition.^ For the two planar 

class, V3xV^ 



classical \/3x \/3 and q—0 GSs we get m ' 



m 



class, g—O 



0.40825 (note that for the kagome 
lattice one obtains the same value). For the quantum 
model one finds (m+)2 = 0.149113 {N = 18); 0.114822 
(24); 0.094831 (30); 0.082299 (36A); 0.079351 (36B); 
0.073251 (42). For comparison we quote the values for 
the iV = 36 kagome lattice: {m+)^ = 0.059128, and the 
iV = 36 triangular lattice: {m+f = 0.124802. 

We have performed finite-size extrapolations based on 
the standard formulas for 2D spin-ha lf Heiscnberg anti- 
ferromagnets (see, e.g. Refs. I7l24l25j) . namely eo(£) = 
Ao + # +0(L-4) for the GS energy per bond, m+{L) = 



Bo 
Go 



L 

G2 



OiL- 
-0{L 



^) for the order parameter, and A(L) = 
"^) for the spin gap, where Aq = eo(oo). 



TABLE II: Results of the finite-size extrapolation of the GS 
energy per bond eo and the order parameter m"*" (eq. Q) 
of the spin-half HAFM for the star [N = 18,24,30,36,42), 
the kagome (A = 12, 18, 24, 30, 36) and the triangular (A = 
24, 30, 36) lattices. To see the effect of quantum fluctuations 
we scale m"*" by its corresponding classical value m\ 



lattice 



eo 



triangular 



-0.1842 
0.386 



kagome 



-0.2172 
0.000 



star 



-0.3091 
0.122 





FIG. 5: Finite-size behavior of the spin gap A, i.e. the gap 
to s = 1 excitations versus inverse system size 1/A''. Results 
are shown for the star lattice (triangles) in comparison with 
those for the kagome lattice (circles, compare Ref. 1231 and 
the triangular lattice (squares). 



Bq = r7i+(oo), Go = A(oo) and L = N^/^. We present 
the results of the extrapolation for eo and 77i+ in compar- 
ison with the triangular and the kagome lattice in Table 
nil The HAFM on the star lattice has lowest GS energy 
eo; the extrapolated order parameter is finite but very 
small. 

Fig. shows the finite-size behavior of the spin gap 
of the star lattice in comparison to the kagome and tri- 
angular lattices. We extrapolate a quite big spin gap 
A = 0.380 for the star lattice. For the triangular lat- 
tice, the data with A^ < 36 seem to suggest a finite spin 
gap which is, however, spurious. This illustrates that the 
extrapolation of the spin gap may be most affected by 
systematic errors.^ Nevertheless, the spin gap of the star 
lattice in Fig. |S1 exhibits only comparably small finite- 
size effects. Hence, the estimation of a non-zero spin 
gap for the star lattice seems to be reliable. The spin 
gap extrapolated for the kagome lattice^ is more than 
six times smaller, but note that the existence of a gap 
for the kagome lattice at all is not fully clear, as is also 
evident from Fig. O 



FIG. 6: Magnetization curves for star lattices with A'^ — 18, 
36A and 36B (complete) and N = 42, 54 and 72 (partial). 
Inset: High-field part of the magnetization curves for N = 
36A, 36B, 54 and 72 sites. 



where a represents a pair of sites i,j corresponding to 
a Jd bond. Thus the product in Q is effectively taken 
over all Jd bonds of the star lattice. In Q the vec- 
tors \(f>'j.{0a)) are spin up (down) states at site i with a 
quantization axis corresponding to the classical planar 
GS considered, i.e. [sm{9^)S^ + cos{9a)S^]\(l)±{0a)) = 
±^\(j)'j.{9a)), where the angles 9a correspond to the re- 
spective classical \/3xV3 or q—0 state. |$) depends 
on the variational parameter t and interpolates between 
a rotationally invariant dimer product state for t — 1 
and an uncorrelated planar \/3xy3 or q—0 state for 
t — 0. Optimizing (^|i/|^) with respect to t we get 
E^^'/hond = ~{Jl + JdJt + Jt)/'^'^Jt- For the sublat- 
tice magnetization we obtain 



M = (*| cos(6'„)S'f -I- sin(6l„)S'f I*) 



V-W 



72 



2Jt 



(4) 



for Jd < Jt- M vanishes with a mean- field critical ex- 
ponent at the symmetric point Jd = Jt- Since such an 
approach tends to overestimate the region of the semi- 
classically ordered state^S^ we may interpret the above 
result as a further indication of a dimerized GS. 



V. VARIATIONAL MEAN FIELD APPROACH 

We discuss briefly a variational approach which was 
successfully applied to describe a quantum phase transi- 
tion between Neel phases and a dimer phaseiSSiSl Let 
us consider a model with different NN bonds Jt and Jd ■ 
The GS shall be approximately described by a variational 
wave function 



I*) 



n 



\^X{9a))\q^_{9a))-t\c^U0c.m'^i9^)) 



(3) 



VI. MAGNETIZATION PROCESS 

Finally, Fig. shows magnetization curves of several 
finite star lattices, where the magnetization m is defined 
as m — 2{S^)/N. Due to computational limitations, only 
the high-field region can be studied for TV > 36. For ex- 
ample for A^ = 42, reliable data are available only for 
iTi > 1/3 (and of course m = 0, see section irv|) . Fur- 
thermore, the lowest parts of the curves for A = 54 
(17/27 < 771 < 19/27) and A^ = 72 (3/4 < 777 < 5/6) 
are based on assumptions concerning the symmetry of 
the GS. 



First, one observes a pronounced zero-field plateau in 
Fig- El corresponding to the spin gap discussed in section 
IIVI Candidates for further plateaux emerge e.g. a.t m = 
1/3, 7/9 and 8/9. The finite-size effects at the boundaries 
of these candidate plateaux are comparably weak for m — 
7/9 (see inset of Fig. EJ such that the numerical evidence 
in favor of a plateau with tti 7^ is strongest in this case. 

The presence of a plateau at m = 1/3 is also plau- 
sible since the star lattice consists of triangles. More 
precisely, Ising-like anisotropics can be argued to give 
rise to a plateau at to = 1/3 due to up-up-down con- 
figurations on the triangles. The number A/'conf. of such 
Ising configurations can be determined by explicit enu- 
meration, yielding e.g. A/'conf. = 123 528, 3 508 392, and 
531 606 684 for the N =42, 54, and 72 lattices, respec- 
tively. This number grows asymptotically approximately 
as A/'conf. oc (1.322 . . .)^. The number of Ising configu- 
rations is much bigger than the corresponding number 
on the kagome lattice at ttt, = 1/3 (see Rcf. 28 and ref- 
erences therein), indicating that the tendency towards a 
disordered GS at to = 1/3 may be stronger on the star 
lattice than on the kagome lattice. 

Just below saturation, we see a jump in the magne- 
tization curve (see inset of Fig. Ej). Indeed, the pres- 
ence of this jump follows from a general construction 
of independent localized magnons for strongly frustrated 
lattices, ^'^^ which in the case of the star lattice live on the 
dodecagons. The expected height 6m = 1/9 of the jump 
for sufhciently large iV is confirmed for iV = 54 and 72 
(see inset of Fig. |S| . Note that the existence of localized 
magnon states also leads to a finite residual T — en- 
tropy at the saturation field h = 5 J/3i2£ and may favor 
a tendency towards a spin-Peierls deformationi^ On gen- 
eral grounds one expects a plateau just below the jump, 
i.e. at the candidate value m — 8/9 mentioned beforCf^ 
although the available data do not allow unambiguous 
confirmation of this plateau. 



VII. DISCUSSION AND CONCLUSION 

Similar as for the kagome lattice the results reported in 
this paper yield indications for a quantum paramagnetic 
GS for the star lattice, too. Whereas this statement is 
well-known for the kagome lattice, the star lattice repre- 
sents a new example for a quantum HAFM on a uniform 
2D lattice without semi-classical GS ordering. However, 
we emphasize that despite the fact that on the classi- 
cal and semi-classical level (LSWT) we have very similar 
physics as for the kagome lattice (i.e. one zero mode, the 
classical GS degeneracy is not lifted) the quantum para- 
magnetic GS for the star lattice is of different nature 



than that for the kagome lattice. The quantum GS for 
the star lattice is characterized by extremely strong NN 
correlation on the dimer bonds (more than 60% larger 
than the NN correlation of the honeycomb lattice having 
the same coordination number z = 3) and a weak NN 
correlation on the triangular bonds (only about 30% of 
the NN dimer correlation and significantly less than the 
triangular NN correlation of the kagome and the triangu- 
lar lattices). The singlet-triplet spin gap is particularly 
large (about six times larger than that for the kagome 
lattice). Although the classical GS exhibits a huge non- 
trivial degeneracy remarkably one does not find low- lying 
singlets within this large spin gap, rather the first sin- 
glet excitation is well above the lowest triplet state. The 
low-lying spectrum as a whole resembles the spectrum of 
weakly coupled dimers.^ All these features support the 
conclusion that the quantum GS of the HAFM on the 
star lattice is dominated by local singlet pairing. This 
dimerized GS represents a so-called explicit valence-bond 
crystal stated which respects the lattice symmetry. 

Although we could expect a gapped quantum paramag- 
netic explicit valence-bond crystal GS for the star lattice 
in case of strong dimer bonds Jd 3> Jt, this should be 
contrasted with other models where an explicit valence- 
bond crystal GS arises in the limit of strong dimer bonds. 
For example, in the simple s = 1/2 Heisenberg bilayer 
model, the picture of weakly coupled dimers is qualita- 
tively correct only for an interlayer exchange coupling 
J± significantly larger than the intralayer coupling J 
(Jj_ > 2. 5 J, see Refs.j26;32 and references therein). Be- 
cause the bilayer Heisenberg model is not frustrated, the 
classical GS does not exhibit any non-trivial degener- 
acy and the GS remains the semi-classical Neel state for 
J± Ki J and all values of si^ By contrast, the quan- 
tum paramagnetic GS appears in the s = 1/2 star lattice 
even in the uniform case Jd — Jt- This difference can 
be attributed to the strong frustration present in the star 
lattice. 

The magnetization curve of the s = 1/2 HAFM on the 
star lattice shows a jump just below saturation and sev- 
eral candidates for plateaux e.g. at to = 1/3, 7/9 and 
8/9 as some typical features of strongly frustrated quan- 
tum spin lattices. Furthermore, low-energy excitations 
present for certain magnetic fields promise large magne- 
tocaloric effects 1^ 
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